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Evidence for Quantized Displacement in Macroscopic Nanomechanical Oscillators 
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We report the observation of discrete displacement of nanomechanical oscillators with gigahertz- 
range resonance frequencies at millikelvin temperatures. The oscillators are nanomachined single- 
crystal structures of silicon, designed to provide two distinct sets of coupled elements with very low 
and very high frequencies. With this novel design, femtometer-level displacement of the frequency- 
determining element is amplified into collective motion of the entire micron-sized structure. The 
observed discrete response possibly results from energy quantization at the onset of the quantum 
regime in these macroscopic nanomechanical oscillators. 
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The quantum mechanical harmonic oscillator is a fun- 
damental example of textbook quantum mechanics 0. 
Its direct experimental realization in truly macroscopic 
mechanical systems is of interest to a wide range of fields 
0- El> Q j which include quantum measurement [3, El j 
quantum computation |l2j, atomic and quantum optics 
[l3| , condensed matter physics 0, 0, E| and gravita- 
tional wave detection [7|, ja l ■ However, despite intense ex- 
perimental efforts l- d Il 7 |. Il8 [ . particularly with nanome- 
chanical structures |9j, lll| , signatures of quantum be- 
havior such as quantum transitions in a macroscopic me- 
chanical oscillator are yet to be observed 0, • 

The essential problem in achieving quantized behav- 
ior in mechanical structures |19| has been the access to 
the quantum regime. Two characteristic time scales, de- 
coherence time and dissipation time, define quantum- 
to-classical crossover. Although decoherence imposes a 
much stricter condition, a necessary requirement for ob- 
serving quantum behavior is given by dissipation or en- 
ergy relaxation (1/Q, inverse quality factor): for a sys- 
tem with Q = 1, the quantum of oscillator energy hf 
is larger than or comparable to ksT. Realization of this 
criterion requires both millikelvin temperatures and giga- 
hertz range frequencies. For example, a nanomechanical 
beam with a resonance frequency of 1 GHz will enter the 
quantum regime at T = 48 mK. For a doubly-clamped 
beam, the fundamental frequency scales as ^E/p(t/L 2 ), 
where E is the Young's modulus, p is the mass density, 
t and L are thickness and length respectively. In typi- 
cal materials like silicon, all dimensions must be in the 
sub- micron range to achieve gigahertz resonance frequen- 
cies. However, if structure dimensions are reduced to in- 
crease the resonance frequency, it naturally increases the 
spring constant k. As the structure becomes stiffer, the 
displacement on resonance, x = FQ/k, decreases for a 
given amplitude of force F. For a gigahertz-range beam, 
the typical displacement is on the order of a femtometer. 
Detecting femtometer displacements is further impeded 
because the quality factor is known to decrease with de- 
creasing system dimensions 



Propelled by the recent advances in nanomcchanics, 
there have been numerous attempts to approach the 
quantum regime in nanomechanical oscillators |l7l lisj 
with low thermal occupation number Nth = ksT/hf. 
The central thrust of this effort has been the devel- 
opment of ultra-sensitive displacement detection tech- 
niques. These include the coupling of the nanomechan- 
ical beam to a single-electron transistor (SET) sensor 
[H illH, a cooper-pair box device HcJ, SQUID sen- 
sor 23|, and piezoelectric sensor [28| as well as optical 
interferometric techniques [2^. Recently, in a 116-MHz 
beam oscillator measured down to a cryostat tempera- 
ture of 30 mK, Knobel et al [2lj have reported displace- 
ment sensitivity one hundred times the standard quan- 
tum limit, which is the limiting factor despite the very 
low occupation numbers they achieve. At a lower oscilla- 
tor frequency of 20 MHz measured down to the oscillator 
temperature of 56 mK, LaHaye et al |31| have achieved 
greater resolution, 4.3 times the standard quantum limit, 
at the cost of larger Nth = 58. Our approach is inherently 
different in that instead of attempting to further improve 
the detection sensitivity, we focus on a novel design of 
the nanomechanical structures, which display gigahertz 
range frequencies with a corresponding displacement in 
the picometer range. 

In this Letter, we report the first observation of dis- 
crete displacement of possible quantum origin in a set 
of nanomechanical oscillators, which resonate at a fre- 
quency as high as 1.49 GHz. At a measured temperature 
of 110 mK, which corresponds to a thermal occupation 
of Nth — * 1) the oscillators demonstrate transitions be- 
tween two discrete positions, with consistent amplitude 
in both magnetic field and time domain sweeps. We ar- 
gue that the wave functions of the two low-lying energy 
levels of the oscillator at Nth 1 result in the observed 
quantized displacement. Furthermore, the nanomechan- 
ical structure truly represents a macroscopic quantum 
system as the quantized displacement involves roughly 
50 billion silicon atoms. 

Our device is an antenna-like structure designed to 
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FIG. 1: a) SEM micrograph of the suspended antenna oscilla- 
tor. It consists of a central Si beam, 10.7 fim long and 400 nm 
wide, as well as two arrays of 500 nm long and 200 nm wide 
paddles on both sides. The total thickness of the structure is 
245 nm, comprised of the device layer of silicon (185 nm) and 
the thermally evaporated gold electrode (60 nm) , colorized in 
yellow, b) Modal simulation of the antenna structure, show- 
ing the low frequency (/ ~ 10 MHz) fundamental resonance 
mode, c) In the high order collective mode, the paddles vi- 
brate at their own natural frequency (/ > 1 GHz), and the 
induced in-phase strain drives the central beam, which acts 
as an amplifier of the paddle motion. 



have coupled but distinct components, pictured in Fig. 1. 
Two rows of small paddles on both sides of a doubly- 
clamped beam serve as frequency-determining elements, 
resonating at the same gigahertz natural frequency. 
When all paddles vibrate in phase, shown in Fig. lc, the 
central beam couples to the paddle motion at the same 
gigahertz frequency and effectively amplifies it, due to its 
large size. 

We have fabricated four identical antenna structures 
from single crystal silicon with e-beam lithography and 
nanomachining. When placed at the center of a 16 
Tesla superconducting magnet in a dilution fridge, each 
structure is driven magnetomotively by a Lorentz force 
Fdr(u) = I(ijj)LB with current I(oj) through the elec- 
trode of length L in perpendicular magnetic field B. The 
resulting motion is given by the harmonic oscillator re- 
sponse function x(u>) — Fd r {^) /[(^o — w + iojQto / Q)m\. 
The beam displacement x in magnetic field B induces 
a voltage V em f(uj) = i£LBujox(oj) in the gold electrode. 
Here £ is a mode-dependent integration constant [3(| • We 
measure V em f using an RF network analyzer. Combin- 
ing the expressions for i^ r (w), x(uj) and the harmonic 
oscillator equation, the magnetomotive expression is: 



V em f(u) 
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A typical low order resonance mode of the antenna 
structure is measured at 21 MHz, as shown in Fig. 2a. 



FIG. 2: Low order resonance mode of the antenna structure 
at 21 MHz with Q = 11000. a) The induced voltage V em f is an 
inverted Lorentzian peak on top of a noise background. Inver- 
sion results from extra phase in the coaxial cable impedance. 
Three representative sweeps show the resonance peak at 5, 6, 
and 8 Tesla. b) As expected from the magnetomotive scheme, 
Eq. 1, the peak height varies as B 2 for constant I& T . c) Plot 
of the linear power dependence of the resonance peak at 2.5 
Tesla field, on a log-log scale. The corresponding displace- 
ment x is linear with driving force Fd r for forces below 1.3 
pN. Nonlinearity ensues for larger Fd r - The effective spring 
constant k e ff = rmoQ measured on resonance lo — > loo, is 330 
N/m. 



With the oscillator at the equilibrium temperature of 60 
mK, the thermal occupation is N t h — 60, corresponding 
to the classical regime. We verify the expected B 2 depen- 
dence of the response on the magnetic field B in Fig. 2b, 
which follows from the magnetomotive relation in Eq.l. 
For small forces, the center-beam displacement x on res- 
onance varies linearly with the driving force Fd r in accor- 
dance with x(u!q) — QFd r {uo)/k e ff (Fig. 2c). The linear 
fit yields an effective spring constant k e ff = 330 N/m. A 
detailed analysis of all the other complex modes and their 
shapes as well as their dependence on structure geometry 
will be given elsewhere. 

Thermal equilibrium of the oscillator with the mix- 
ing chamber of the dilution cryostat occurs primarily 
through thermal phonon exchange between the single- 
crystal oscillator and the bulk and by coupling to the gold 
electrode on top of the beam through the Kapitza resis- 
tance 34). The thermal phonon wavelength is given by 
Xth = fiVg/kBT, where v s = 5000 m/s is the sound veloc- 
ity in silicon. Free thermal phonon propagation through 
the silicon beam (width w — 400 nm) can occur down to 
T = hv s /ksw = 95 mK. Thermalization of the oscillator 
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FIG. 3: Classical response of the antenna structure in the 
1.48 GHz collective high order mode with Q — 150. At 
T = 1000 mK the occupation factor is Nth — 14 for this mode. 

a) Four plots (A,B,C, and D) of the induced voltage V em j 
show the main peak at 1.479 GHz and a smaller secondary 
peak at 1.49 GHz, for various values of the magnetic field B. 
The data is normalized to the noise background at Tesla. 

b) In agreement with the magnetomotive expression, the plot 
of the peak height versus magnetic field B shows quadratic 
dependence, entirely analogous to the 21 MHz plot, c) Lin- 
ear dependence of the induced power Pi n on driving power 
Pdr shows direct evidence of classical Hooke's law oscillator 
behavior. The effective spring constant is k e ff = 188 N/m. 



occurs also through the silicon-gold interface. The ther- 
mal Johnson noise signal size St h — \ 4k B T Q decreases 

dramatically with increasing frequency. For example, a 
l-/xm long silicon beam with a resonance frequency of 1.5 
GHz has Si 11 ~ 3 x 10~ 18 m/VHz at a temperature of 
100 mK, three orders of magnitude smaller than the best 
sensitivity of ~ 10~ 15 m/V Hz reported to date |3(X l3l|. 

We observe a high-frequency collective mode at 1.48 
GHz with Q ~ 150, as shown in Fig. 3a; this is the 
highest nanomechanical resonance frequency reported to 
date |33j. For temperatures below 100 mK, the occu- 
pation factor is N t h — > 1, and we expect the onset of 
non-classical behavior. However, it is imperative to ver- 
ify that the mode is classical for larger values of Nth- 
Magnetic field and power sweeps T — 1000 mK, which 
corresponds to Nth — 14, indeed demonstrate the same 
B 2 field and linear force dependence, as shown Figs. 3b 
and 3c. Furthermore, the experimentally determined ef- 
fective spring constant fc e // = 188 N/m is an order of 
magnitude lower than the estimated value for a 1-GHz 
straight silicon beam, giving higher displacements as ex- 
pected from our structure design. 
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FIG. 4: Non-classical response in the 1.49 GHz collective 
mode at T — 110 mK. a) Four discrete plots, showing rapid 
peak growth within the narrow B field range from 6.0 to 6.87 
Tesla, which corresponds to 45 to 52 pN of force. We observe a 
0.6 percent frequency blue shift with decreasing temperature, 
b) The peak growth is non-monotonic, as seen in this separate 
continuous field sweep at the center frequency fo = 1.49 GHz. 
The peak undergoes several transitions with amplitude 500 
nV between two well-defined states. The labels A,B,C, and 
D indicate the approximate location of peak growth in plot 
a), c) Two sweeps of the magnetic field, up (black) and down 
(blue), reproducing the same transition, with occasional tran- 
sitions to a possible intermediate state, d) With all parame- 
ters held constant, the system spontaneously decays in time, 
with no further changes observed within a measurement pe- 
riod of 2 hours. 



In stark contrast to the classical behavior at T = 
1000 mK, the collective mode response exhibits non- 
monotonic magnetic field dependence when the tem- 
perature is lowered to T = 110 mK, corresponding to 
Nth — 1- In Figs. 4b and 4c, the plots show discrete peak 
voltage transitions at three values of the magnetic field. 
While the transitions do not always occur at exactly 
the same field values, the size of the jumps remains un- 
changed, suggesting that the oscillator switches between 
two well-defined states. The transition is reproducible for 
field sweeps in either direction. Furthermore, in prelim- 
inary measurements we have observed spontaneous de- 
cay from the upper to the lower state in time, shown in 
Fig. 4d, with a more detailed investigation forthcoming. 

The oscillator response in the 1.48-GHz mode in Figs. 3 
and 4 demonstrates continuous dependence on drive force 
(energy) at 1000 mK and discrete response at 110 mK, 
corresponding to the classical regime (Nth ~ 14) and 
the quantum regime (Nth ~ 1) respectively. Since our 
experiment involves the measurement of the absolute 
value of the displacement |a;|, it is possible that the dis- 
crete transitions show signatures of energy quantization 
as \x\ ex y/E, where E is the oscillator energy. The two 
distinct saturated states (A and D in Fig. 4b) could then 
correspond to the ground state |0) and the first excited 
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state 1 1) in the energy eigenbasis, respectively. The in- 
termediate jumps between these two states (in the range 
of 6.5-7 tesla in Fig. 4b) could represent transitions in- 
duced by thermal fluctuations, as the temperature even 
at Nth ~ 1 is high enough that the thermal energy fc^T 
smears the gap hf in the energy levels. If this is true, 
then future experiments at higher frequencies and lower 
temperatures, corresponding to a lower Nth, will show 
clear transitions from |0) — » |1). Further reduction in Nth 
could even show multiple transitions |0) — * |1) — * |2) • • •, 
corresponding to access to the first and higer-order ex- 
cited states. 

Another possibility is that the oscillator does not start 
from the energy eigenbasis, as there is no a priori reason 
for it to be in this preferred basis. As the driving energy 
is increased (by increasing the field), the oscillator goes 
from a certain ground state (point A in Fig. 4b) to a 
linear combination of |0) and |1) with coefficients corre- 
sponding to maximum amplitude of displacement. This 
process should necessarily involve only the first couple of 
low-lying states. In fact, for Nth ~ U only the low-lying 
states (|0) and |1)) are substantially populated. 

Similar discrete behavior can result classically due to 
transitions between two nonlinear bistable beam states. 
After an extensive study of bistability in nanobeams |35| . 
we rule out this mechanism in this case because the re- 
sponse of the structure at the drive level used in the ex- 
periment is manifestly linear, seen from the Lorentzian 
peak shape (Fig. 4a). However, there may be other semi- 
classical mechanisms which involve both nonlinearity and 
quantum effects pSfij . 

In order to clearly establish whether or not the ob- 
served discrete features are quantum mechanical in ori- 
gin, a complete theory is needed. For example, simula- 
tion of the mechanics of the antenna structure containing 
fifty billion silicon atoms must involve multiscale model- 
ing, including molecular dynamics at the atomic level 
and finite element analysis in the continuum elastic the- 
ory. Understanding the measurement of the quantum 
system itself requires identifying the components of the 
quantum measurement process: quantum system, mea- 
suring apparatus, and their interaction. Coupling to the 
environement, which introduces decoherence and dissi- 
pation, must be included in the theoretical analysis. On 
the experimental side, an irrefutable observation of quan- 
tum mechanical transitions will require measurements at 
higher frequencies and lower temperatures (Nth < 1) 
with a clear gap between the energy levels. Detailed 
study of the noise in the discrete states will also help 
elucidate the possible quantum mechanical origin of the 
oscillator response. 

In conclusion, we have measured the displacement in 
a series of nanomechanical oscillators with 1.4- GHz reso- 
nance frequencies down to 110 mK temperatures. While 
the low frequency mode at 21 MHz shows classical be- 



havior with expected drive dependence (oc B 2 ), the 1.49 
GHz mode displays non-monotonic dependence on driv- 
ing force at a temperature which corresponds to a ther- 
mal occupation number N t h — > 1. Our experimental data 
indicates the first observation of quantum displacement 
in macroscopic nanomechanical oscillators. Our work is 
supported by the NSF, ARL, ACS and the Sloan Foun- 
dation. 
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